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Linear Stability of Equilibrium Points in the 
Generalized Photogravitational Chermnykh's 
Problem 



Badam Singh Kushvah 1 



OO 

O 

o 

(N 

c 



c/3 
Q 

o3 



> 
(N 



O 

OO 

o 



X 



© Springer- Verlag •••• 

Abstract The equilibrium points and their linear sta- 
bility has been discussed in the generalized photograv- 
itational Chermnykh's problem. The bigger primary 
is being considered as a source of radiation and small 
primary as an oblate spheroid. The effect of radi- 
ation pressure has been discussed numerically. The 
cohinear points are linearly unstable and triangular 
points are stable in the sense of Lyapunov stability pro- 
vided fjb < fXRouth = 0.0385201. The effect of gravita- 
tional potential from the belt is also examined. The 
mathematical properties of this system are different 
from the classical restricted three body problem. 

Keywords Equilibrium Points: Linear Stability: Gen- 
eralized Photogravitational: Chermnykh's Problem: 
Radiation Pressures 



1 Introduction 

The Chermnykh's problem is new kind of restricted 
three body problem which was first time studied by 
Chermnykh (1987). Recently many authors studied 
this problem such as Jiang and Yeh (2004a,b,c) con- 
sidered the influence from the belt for planetary sys- 
tems and found that the probability to have equilib- 
rium points around the inner part of the belt is larger 
than the one near the outerpart. Papadakis (2005) ex- 
amined the motion around the triangular equilibrium 
points of the restricted three-body problem under an- 
gular velocity variation. Yeh and Jiang (2006) studied 
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a Chermnykh-Like problem in which the mass param- 
eter \i is set to be 0.5. Jiang and Yeh (2006) found 
the equilibrium points in the Chermnykh-Line problem 
when an addition gravitational potential from the belt 
is included. The solar radiation pressure force F p is ex- 
actly apposite to the gravitational attraction force F g 
and change with the distance by the same law it is pos- 
sible to consider that the result of action of this force 
will lead to reducing the effective mass of the Sun or 
particle. It is acceptable to speak about a reduced mass 
of the particle as the effect of reducing its mass depends 
on the properties of the particle itself. 

Ishwar and Kushvah (2006) examined the linear sta- 
bility of triangular equilibrium points in the generalized 
photogravitational restricted three body problem with 
Poynting- Robertson drag, L4 and L5 points became un- 
stable due to P-R drag which is very remarkable and 
important, where as they arc linearly stable in classical 
problem when < \i < ^R OU th = 0.0385201. Kushvah, 
Sharma, and Ishwar (2007a, b,c) examined normaliza- 
tion of Hamiltonian they have also studied the nonlin- 
ear stability of triangular equilibrium points in the gen- 
eralized photogravitational restricted three body prob- 
lem with Poynting-Robertson drag, they have found 
that the triangular points are stable in the nonlinear 
sense except three critical mass ratios at which KAM 
theorem fails. Papadakis and Kanavos (2007) given nu- 
merical exploration of Chermnykh's problem, in which 
the equilibrium points and zero velocity curves studied 
numerically also the non-linear stability for the triangu- 
lar Lagrangian points are computed numerically for the 
Earth-Moon and Sun- Jupiter mass distribution when 
the angular velocity varies. The stability of triangle 
libration points in generalized restricted circular three- 
body problem has been studied by Beletsky and Rod- 
nikov (2008). Das, Narang, Mahajan, and Yuasa (2008) 
examined the stability of location of various equilibrium 
points of a passive micron size particle in the field of 
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radiating binary stellar system within the framework of 
circular restricted three body problem. 

In this paper we have obtained the equations of 
motion, the position of equilibrium points and their 
linear stability in the generalized photogravitaional 
Chermnykh's problem. The effect of radiation pres- 
sure, oblatcness, and gravitational potential from the 
belt has been examined analytically and numerically. 
We have seen that the collinear equilibrium points are 
linearly unstable while the triangular points are condi- 
tionally stable. 



2 Equations of Motion and Zero Velocity 
Curves 

We consider the baryecntric rotating co-ordinate sys- 
tem Oxyz relative to incrtial system with angular ve- 
locity u) and common z-axis. We have taken line join- 
ing the primaries as x-axis. Let mi, m 2 be the masses 
of bigger primary(Sun) and smaller primary(Earth) re- 
spectively. Let Ox, Oy in the equatorial plane of 
smaller primary and Oz coinciding with the polar axis 
of m 2 . Let r e , r p be the equatorial and polar radii of m 2 
respectively, r be the distance between primaries. Let 
infinitesimal mass m be placed at the point P(x, y, 0). 
We take units such that sum of the masses and distance 
between primaries is unity, the unit of time i.e. time 
period of mi about m 2 consists of 2ir units such that 
the Gaussian constant of gravitational k 2 = 1. Then 
perturbed mean motion n of the primaries is given by 

2 _ 2 

n 2 = 1 + ^4 a i where A 2 = r \Ji v is oblateness coef- 



ficient of m 2 . Let u 
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with mi > m 2 , where fi is mass parameter. Then co- 
ordinates of mi and m 2 are (— /z,0) and (1 — /z, 0) re- 
spectively. In the above mentioned reference system we 
determine the equations of motion of the infinitesimal 
mass particle in xy-plane as Kushvah (2008). 



x-2ny = U x , 
y + 2nx = U y 



(1) 
(2) 



where 
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where 



= n 2 (x 2 +y 2 ) (1 - /z)gi /z ^2 



gi = 1 — -p£ is a mass reduction factor expressed 
in terms of the particle radius a, density p radiation 
pressure efficiency factor \ (i n C.G.S. system): q\ = 



5.6x10' 
ap 



■X < 1 • The assumption q\ 



constant 

is equivalent to neglecting fluctuations in the beam of 
solar radiation and the effect of the planets shadow. 

2.1 Miyamoto and Nagai (1975) Profile Model 

In this model we introduce the potential of belt as: 



V{r,z) = - 



(4) 



x 2 +y 2 , 



^r 2 + (a+Vi 2 ~+b 2 ") 2 

where Mb is the total mass of the belt and r 2 = 
a, b are parameters which determine the density profile 
of the belt. The parameter a controls the flatness of 
the profile and can be called "flatness parameter" . The 
parameter b controls the size of the core of the density 
profile and can be called "core parameter" . When a = 
b = the potential equals to the one by a points mass. 
In genrel the density distribution corresponding to the 
above V(r, z) in ( 4) is as in Miyamoto and Nagai (1975) 



p(r,z) 

where N - 
obtained 

V(r,0) = 
and V x = 



b 2 M b [ar 2 + (a + 3zV)] (a + N'f 
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r 2 + (a + N) 

a + b, z 
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(5) 



Vr 2 + T 2 



M b x 



v v 



MbV 



0. Then we 
(6) 

Now we con- 



( r 2 +T 2)3/2 ! *y — ( r 2 +T 2j3/2 1 

sider only the orbits on the x — y plane, then the equa- 
tions of motion are modified by using ( 1,2) in the fol- 
lowing form: 



x — 2ny 
y + 2nx 



U x 



v x 



(7) 
(8) 



where 
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Then the perturbed mean motion n of the primaries 
is changed into the form n 2 = 1 + ^j 2 - + ^ 2 ^^ 3 y 2 , 

where r 2 = (1 — fi)qi^ 3 + /J 2 , we set r = r c = 
0.08, T = 0.01, n = 0.025 for further numerical re- 
sults. The energy integral of the problem is given by 
C = 217 — x 2 — y 2 , where the quantity C is the Jacobi's 
constant. The zero velocity curves[see figure ( 1)] are 
given by: 



C = 2Q(x,y) 

3 Position of Equilibrium Points 



(10) 



The position equilibrium points of Chermnykh's pro- 
belm is given by putting Q x = Cl y = i.e., 



n 2 x 



1 - fi)qi(x + n) ii(x + ji-l) 



3 \iA 2 {x + 11—1) 
2 
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M b y 



+ T 2 ) 



3/2 



0,(11) 



(12) 



From equation ( 12) n 2 - _ 4 3M2 



AL, 



or y = 0. 



( r 2 +T 2)3/2 

3.1 Collinear Equilibrium Points 
In this case suppose 

(1 - fi)qi(x + fi) (i(x + fi-l) 



f(x,y) = n 2 x 

1 1 

3fiA 2 (x + /i- 1) M b x 
2 r\ < r 2 + T 2 



3/2 



0,. 



(13) 



g(x,y) 



[1-Li)qi fl Z\lA 2 



(r 2 +T 2 f /2 
f(x,0) = P(x) + Q(x) 



y = 



(14) 
(15) 




Fig. 1. — The zero velocity curves (l):M b = 0, 
(2):M 6 = 0.2,(3):M 6 = 0.4, (4):M b = 0.6, when 
li = 0.025, r = 0.8, T = 0.01, A 2 = 0, q = 1 
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Fig. 2. — Figure show the position of the collinear 
equilibrium points in frame (a) when q\ = 1 , A 2 = 
and frame (b) when qi = 1,A 2 = 0.02, the curves 
are leveled by (1-4) correspond to the mass of belt 
M b = 0,0.2,0.4,0.6. 
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where 



(2008) then from equations ( 7) and ( 8) we obtained: 



n 2 x — 



P(x) 
3 fiA 2 {x + /i 



(1 - n)qi(x + (j) n(x + fj,-l) 



\x + 



1) 



2 |z + /i-l| 5 



Q(x) 



M b x 



(x 2 + T 2 ) 
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\x + fi — 1| 3 
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To investigate the position of collincar equilibrium 
points divide the orbital plane Oxy into three parts 
with respect to the primaries x < —fi, 1 — \i < x and 
— fi < x < 1 — (i, for each part the function P(x) is 
defind as follows: 



3 {rl + T 2 ) 3/2 
From above, the triangular cquilibriumpoints are as: 



(19) 
(20) 



x = -jj, + 
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fJ,A 2 
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(z+m) 2 

I 3 M2 
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If X < — (i, 



If 1 — jj, < X, 



If — [J, < X < 1 — fj, 

(18) 



+- 



M (2r c - l)M fc A 2 



(1- M ) (r c 2 +T 2 ) 



3/2 



(21) 



2/3 



(4 - , 2 / 3 ) + 2 (, 2 / 3 - 2) A 2 



4(2r c - l)M t 



(^ /3 - 3) 



3nA 2 (ql /3 -3) 
2(1-M) 



3(r c 2 +T 2 ) 



3/2 



-,1/2 
22) 



When a; G (-co, -/z), lim^-oo /(or, y) < 0, 
lim 2 ,^_ AI - f(x,y) > 0, f"(x) > 0, Q'(x) < so there 
is a point xs for which f(x3, 0) = 0. If x € (1 — /i, oo), 
lim x ^ (1 _ /j)+ /(a;, j/) < 0, ]im x ^ 00 f(x,y) > 0, P'{x) > 
0, Q'(x) > so there is a point x 2 for which /(x 2 , 0) = 

0. Now x e (— 1 - A*): lim x ^ + f(x, y) > 0, 
lim^n^)- f{x,y) > this implies that an even(or 
zero) number real roots of f(x, y) exists in this range. If 
x e (— /i, 0), consider two cases (i) x € (— (i, — and 
(ii)xe (-^), 0),we obtained Q{-T_)+p{-T_) > 

1. e. i(x, 0) 7^ so there is no equilibrium point at 
x = If Q(-^) + i > (-^) > 0, we obtained 
lim x ^ - f(x,y) < 0, \\m x ^_^+ f{x,y) < so there 
exists two new equilibrium points x\,\ G ( — ^tjjO) and 

■x b2 e (-M, -^) for which f(x, 0) = 0. ButQ(-^) + 
P{—-^) < the there exists no equilibrium points in 

(-M,0). IfT< V2 M andQ(-^) +P(-^) > then 
we have two new equilibrium points. Hence we have 
found there are five equilibrium points on the £-axis 
for the given system. The position of above points are 
presented graphically by frames a, b in figure ( 2) the 
curves are leveled by (1-4) correspond to the mass of 
belt M b = 0,0.2,0.4,0.6. 

3.2 Triangular Equilibrium Points 

The triangular equilibrium points are given by putting 
£l x = £l y = 0, y 0. Using the method as Kushvah 



All these results are similar with Szebehely (1967), Ra- 
gos and Zafiropoulos (1995), Jiang and Yeh (2006) and 
others. 



4 Linear Stability 

To study the linear stability of any equilibrium point 
change the origin of the coordinate system to its posi- 
tion (x*,y*) by means of x = x * +a, y = y * +/3, 
where a = £e At , j3 = rje xt are the small displacements 
£, rj, A these parameters, have to be determined. There- 
fore the equations of perturbed motion corresponding 
to the system of equations ( 7), ( 8) may be written as 
follows: 



a - 2n/3 = afl* xx + f3tt* xy 
j3 + 2na = afi; + f3tt* yy 



(23) 
(24) 



where superfix * is corresponding to the equilibrium 
points. 








(\ 2 -n xx )z + (-2n\-n xy ) v 
(2n\-n; x )t + (\ 2 -n; y )r, 

Now above system has singular solution if, 
A 2 - ft* 



(25) 
(26) 



XX 

2n\ - Q 



-2n\-Q xy 
A 2 - ft* 

yx n L yy 
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A 4 + 6A 2 + d = 



(27) 



At the equilibrium points equations ( 1),( 2) gives us 
the following: 



b = 2n 2 - f, 



3nA 2 



r: 



+ 



d = 9 M (1 - »)yl 

3M b 



3M b T 2 

ir 2 +T 2)5/2 

<ll 



' 1*' 2* 



rp2 



.5/2 



+ 



(l-M)(l + £ft 



3 A 2 
271 



+ 



where /, = &=JpSL + -g- (l , ^ , , „ 

The points L\,L2,Lz, xn, x b2 lie along the line join- 
ing the primaries, from ( 19) we have r\ 
case 



3M b 



1/3 

In this 



(1 - n)q 1 



+ 



x + A 1 ! 3 
3M b 



+ 



+ M -1| 3 



1 + 



. J"2 



5/2 



3 A 2 
2|x + A <-l| 2 

(28) 



With the help of equation ( 28) we have drawn fig- 
ure 3, the different curves (l)-(4) correspond to M b = 
0.0, 0.2, 0.4, 0.6. From this figure we observe that for 
each L\,L%,Lzi x b \ and x b 2, f > 1, in this case charac- 
teristic equation ( 27) has at least one positive root, this 
implies that collinear equilibrium points are unstable in 
the sense of Lyapunov. 




Fig. 3. — x — f, curves correspond to (l):M b = 0, 
(2):M b = 0.2,(3):M 6 = 0.4, (4):M 6 = 0.6, when fj, = 
0.025, r = 0.8, T = 0.01, A 2 = and q x = 1 

Now we have to study the linear stability of trian- 
gular equilibrium points, in this regard we obtained 



b = 9/i(l — fJ,)g, where g 



+ 



3M b 



Ml 



+ 



(1 "">( 1+ *) 



tcristic equation ( 27) we obtained: 



i 



A 



-b±(b 2 -Ad) 1 ' 2 



From charac- 



(29) 



For stable motion < Ad < b 2 , i.e. 
(n 2 -3/M 2 ) 2 >36/i(l-/i)g 

In classical case A 2 = 0,M b = 0, gi = 1, n = 1, we 
have following: 1 > 27^(1 - n) => fj, < 0.0385201. 
Using equation ( 29) we obtained imaginary roots 
Ai,2 = iiwi, A3. 4 = ±iw2, i = V~- T. The characteristic 
frequencies 101,2(0 < 1^2 < are presented by frames 
(a)-(d) of ( 4- 7) in parameter plots u>i — qi for different 
values of A 2 ,qi,M b . They are given in table( 1). We 
observe that they are decreasing function of radiation 
pressure and increasing functions of A2 , M b . Hence the 
triangular equilibrium points are stable in the sense of 
Lyapunov stability provided /1 < fiRouth = 0.0385201. 
In this case we obtained the three main cases of 



resonances: 



u>i — ku>2 = 0, 



1,2,3 



(30) 



For k = 1 we have positive stable resonance and for 
k = 2, 3 we have unstable resonances. Using ( 29) and 
( 30) we obtained a root of mass parameter: 



3g + 2Khb 2 - ^gy/9g - 4Kb 2 + 12^2 , 01 , 
= 777 — iTIal ( 31 ) 



6(g + Kb 2 ) 



where K 
b 2 = A 2 



JW+lJ 2 



n 2 + 



2rM b 



3M b T 2 



1 , 5(2r-l)A/ b 
^ (r 2 +T 2 ) 3 / 2 



( r 2 +T 2)3/2 I ( r 2 +T 2)5/2) 

Now we suppose q\ = 1 — e, 

with |e| << 1, neglecting higher order terms, we ob- 
tained the critical mass parameter values corresponding 
to k = 1, 2, 3 as : 

Mi = 0.0385208965 + 0.0375419787^ 2 

-0.0089174706e - 0.0678734040M 6 (32) 

Ai 2 = 0.0242938971 + 0.0254350205A 2 

-0.0055364958e - 0.0421398438M 6 (33) 

pa = 0.0135160160 + 0.0148764140A 2 

-0.0030452832e - 0.0231785159Af fc (34) 

The linear stability region and main resonance curves 
k = 1,2,3 are shown by /1 — q\ parameter space frames 
(a)-(d) in figures ( 8, 9). The curve corresponding to 
k = l,(qi = l,A 2 = 0,M b = 0,/ii = URouth = 
0.038521) is actual boundary of the stability region. 



D 2 0.4 6 8 1 0.2 4 0.6 0.8 1 p 02 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

V 11 » 11 




Fig. 4. — loi — qi, frames correspond to (a):M b = 0, 
(b):M b = 0.2,(c):Af b = 0.4, (d):M b = 0.6, curves (1) : 
A 2 = 0, (2) : A 2 = 0.02, (3) : A 2 = 0.04, (4) : A 2 = 
0.06, ^ = 0.025, r = 0.8, T = 0.01 



Fig. 7.— w 2 - ?i, frames (a):M 6 = 0, (b):M b 



0.2,(c):M b = 0.4, (d):M b 
0.8, T = 0.01 



0.6, when /i = 0.025, r 




Fig. 5. — 0J2 — 0i, frames correspond to (a):M b = 0, 
(b):M b = 0.2,(c):M b = 0.4, (d):M b = 0.6, curves (1) : 
A 2 = 0, (2) : A 2 = 0.02, (3) : A 2 = 0.04, (4) : A 2 = 
0.06, n = 0.025, r = 0.8, T = 0.01 
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Fig. 8. — (J, — 0i frames correspond to (a):M b = 
(b):M b = 0.2,(c):A/ b = 0.4, (d):M b = 0.6, when A 2 = 
0,k = 1-3, r = 0.8, T = 0.01 
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Fig. 6.— wi - 0i frames (a):M b = 0, (b):M b = 
0.2,(c):M b = 0.4, (d):M b = 0.6, when // = 0.025, r = 
0.8, T = 0.01 



Fig. 9 

(b):M b 
0.02, k 



. — (J, — q\ frames correspond to (a):M b = 
= 0.2,(c):A/ 6 = 0.4, (d):M b = 0.6, when A 2 = 
= 1-3, r = 0.8, T = 0.01 
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The critical values of mass parameter fik are presented 
in table ( 2) for various values of qi,A2,Mb- The classi- 
cal critical values of fi are similar to Dcprit and Dcprit- 
Bartholome (1967). These results are similar to the 
results of Markcllos, Papadakis, and Perdios (1996), 
Kushvah (2008) and others. We observe that the ef- 
fect of radiation pressure reduces the linear stability 
zones and the fik is an increasing function of Mb, A 2 . 

5 Conclusion 

The points L\, L 2 , L3, Xbi, Xb2 he along the line joining 
the primaries. We observe that the effect of radiation 
pressure reduces the linear stability zones, these are also 
affected by belt and the oblateness of second primary. 
The collinear equilibrium points arc unstable while tri- 
angular equilibrium points are stable in the sense of 
Lyapunov stability provided \x < fiR OU th = 0.0385201. 
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Table 1 wi, 2 when r = 0.8, T = 0.01, n = 0.025 



A 2 


gi 


wi : M b = 


u> 2 : M b = 


u)i : M b = 0.2 


u 2 : M b - 0.2 


o)i : M b = 0.4 


U2 : M b = 0.4 


o)i : Mb = 0.6 


uj 2 : Mb = 0.6 


0.0 


1.0 


0.890141 


0.455686 


1.18033 


481 5937 


1.41795 


489382 


1.62232 


4931 27 




0.75 


0.880622 


0.47382 


1.17804 


0.487083 


1.41737 


0.491041 


1.62238 


0.492928 




0.5 


0.869076 


0.494679 


1.17602 


0.491953 


1.41733 


0.491159 


1.62304 


0.490779 




0.25 


0.853749 


0.520684 


1.17436 


0.495895 


1.41812 


0.488893 


1.62461 


0.485532 




0.0 


0.821584 


0.570088 


1.17349 


0.497955 


1.4215 


0.478958 


1.62926 


0.4697 


0.02 


1.0 


91 0983 




1.1934 


477057 


1.42768 




1 63005 


0.491211 




0.75 


0.901845 


0.463871 


1.19127 


0.482336 


1.42716 


0.48813 


1.63013 


0.490934 




0.5 


0.891743 


0.483005 


1.18941 


0.486917 


1.42716 


0.488124 


1.6308 


0.488709 




0.25 


0.878605 


0.506511 


1.18791 


0.490552 


1.42798 


0.485737 


1.63238 


0.48339 




0.0 


0.852388 


0.549485 


1.18724 


1.43137 


1.43137 


0.475658 


1.63701 


0.467476 


0.04 


1.0 


0.929538 


0.439272 


1.20624 


0.472778 


1.43732 


0.483883 


1.63772 


0.489328 




0.75 


0.921985 


0.45491 


1.20426 


0.477788 


1.43685 


0.48529 


1.63783 


0.488972 




0.5 


0.913034 


0.47262 


1.20254 


0.482095 


1.43689 


0.485164 


1.63851 


0.486674 




0.25 


0.901559 


0.494157 


1.2012 


0.485439 


1.43774 


0.482658 


1.6401 


0.481283 




0.0 


0.879387 


0.532615 


1.2007 


0.486672 


1.44113 


0.472436 


1.64471 


0.465288 



Note. — Table 1 presents the roots of characteristic equation( 27). 



Table 2 fi k {A 2 ,M h ), when r = 0.8, T = 0.01 



91 


k 


Hk (0,0) 


)i k (0,0.02) 


W (0,0.4) 


Hk (0,0.6) 


M0. 02,0) 


Mfc(0.02,0.2) 


MO.02,0.4) 


MO.02,0.6) 


1.0 


1 


0.0385209 


0.0525812 


0.0688051 


0.0861218 


0.0404877 


0.0539744 


0.0696964 


0.0863953 




2 


0.0242939 


0.0329695 


0.0428408 


0.0532015 


0.0255597 


0.0339343 


0.0435889 


0.0537117 




3 


0.013516 


0.0182676 


0.0236236 


0.0291855 


0.0142309 


0.0188392 


0.0241121 


0.0295953 




4 


0.00827037 


0.0111565 


0.014396 


0.0177443 


0.00871096 


0.0115163 


0.0147154 


0.01803 




5 


0.0055092 


0.00742441 


0.00956953 


0.0117815 


0.0058038 


0.00766762 


0.00978937 


0.0119837 


0.75 


1 


0.0363201 


0.051579 


0.0684542 


0.0863238 


0.0382382 


0.0529886 


0.0693753 


0.0866222 




2 


0.0229262 


0.0323547 


0.0426289 


0.0533212 


0.024159 


0.0333263 


0.0433931 


0.0538481 




3 


0.0127632 


0.0179323 


0.0235093 


0.0292495 


0.0134588 


0.0113141 


0.0240058 


0.0296688 




4 


0.0078121 


0.0109532 


0.014327 


0.0177827 


0.00824062 


0.0113141 


0.014651 


0.0180743 




5 


0.00520474 


0.00728962 


0.00952388 


0.0118069 


0.00549121 


0.0075334 


0.00974673 


0.012013 


0.5 


1 


0.0341355 


0.0507482 


0.0685365 


0.0872546 


0.0359977 


0.0521785 


0.0694903 


0.0875708 




2 


0.0215661 


0.0318447 


0.0426786 


0.0538727 


0.0227616 


0.0328263 


0.0434633 


0.054418 




3 


0.0120136 


0.0176539 


0.0235361 


0.0295437 


0.0126877 


0.0182322 


0.0240439 


0.0299757 




4 


0.00735548 


0.0107843 


0.0143431 


0.0179594 


0.00777057 


0.0111476 


0.0146741 


0.0182594 




5 


0.00490128 


0.00717768 


0.00953459 


0.0119234 


0.00517872 


0.00742292 


0.00976201 


0.0121354 



Note. - 



Tabic 2 is presents the roots of characteristic equation( 27) 



